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Confinement induced resonances (CIR) in quasi-lD systems have been theoretically predicted and 
observed in various ultracold atomic gases. Here a regularized local frame transformation method is 
developed to treat CIR in a quasi- ID system that has an arbitrary transverse trap shape in general. 
The method is applied to predict the CIR position in a system whose transverse geometry is a square 
well. 



o 

U 

Oh 
< 

O 



43 

Oh 
i 



o 

Oh 



> 

o 
oo 

O 
m 

> 



I. INTRODUCTION 

Quantum systems in reduced dimensions have at- 
tracted extensive interest dating back at least to the 
1950s. Since then, striking predictions have been made 
for the behavior of quantum gases in reduced dimensions. 
In both many-body [TH5] and few-body physics [BHH] , pre- 
dictions exist that quantum states in reduced dimen- 
sional systems are predicted to show qualitatively dif- 
ferent behavior from isotropic 3D systems, both in the 
zero-temperature limit and in the thermodynamic limit. 
Many of these conclusions can be connected to theoret- 
ical explorations from the mathematics community into 
the nature of scattering in a space of arbitrary dimen- 
sions [5], dating even farther back in time. However, 
for a long time, the creation and observation of reduced 
dimensional systems were hampered by limits on exist- 
ing experimental techniques. Later, after laser cooling 
and trapping techniques of ultracold atomic gases were 
developed [10], many experimental quasi-lD and quasi- 
2D [TTHLljj quantum systems were realized and a number 
of theoretically-predicted novel behaviors were observed. 
In quasi-lD fermionic 40 K [T^HS] and bosonic 87 Rb [12] . 
various properties from pure-ID theory have been con- 
firmed. By tuning the ratio(s) of trapping frequencies in 
different directions {lo p /lo z e.g.), the shape of an ultra- 
cold gas can be engineered and varied continuously from 
a nearly isotropic 3D geometry to a quasi-lD or quasi-2D 
geometry. Cross-dimensional effects have been observed 
QUOl] and also theoretically studied [SJQ!]. These phe- 
nomena in low-dimensional systems have recently been 
reviewed by [20H24]. 

After achieving success in creating reduced- 
dimensional systems, physicists began to explore 
the scattering properties of atoms, hoping to find addi- 
tional tunability of interactions in ultracold gases. The 
confinement-induced resonance(CIR) is one such useful 
tunability of particular interest. From an experimental 
point of view, the discovery of CIR adds another way 
to manipulate and control low-dimensional systems. 
Dunjko et.ai [2T] has reviewed developments of CIR for 
a variety of systems. The CIR is unique in the sense that 
the system shows unitary behavior at a non-resonant 
value of the 3D scattering parameters (i.e. the scattering 
length for s-wave, and the scattering volume for p-wave). 



The CIR in a harmonically transverse-confined quasi- ID 
system was predicted theoretically by Olshanii [25], and 
later extended to fermions by Granger and Blume [26] . 
Olshanii [25] derived the exact wave function and phase 
shift for a regularized zero range potential that are asso- 
ciated with low energy s-wave scattering. This problem 
can also be treated using Green's function methods. 
Implementation of an eigenfunction expansion of the 
Green's function for treating a bound state induced by 
a s-wave zero-range potential was discussed in [37]. A 
similar eigenfunction expansion of Green's functions 
is embedded in Olshanii's original work, although not 
explicitly emphasized. Later, a systematic discussion 
of the Green's function method and the corresponding 
Lippman-Schwinger equation was presented in [28] 129], 
Recently the Green's function method was applied 
to treat s-wave zero-range interaction in asymmetric 
transverse harmonic confinement [30] . However for 
higher partial waves, different mathematical models of 
the zero-range interaction are needed in the Olshanii 
treatment [31] . In an alternative treatment, Blume 
[26] used the frame transformation method, in which 
the 3D phase shift information directly determines the 
ID reaction matrix K. Thus the frame transformation 
approach avoids the need to design a zero-range model 
potential for higher partial wave scattering, which has 
some conceptual advantages. For example, a p-wave 
CIR was predicted, later observed experimentally [3"2"] . 
In realistic atomic interactions, the incoming wave in the 
z-direction will not only couple to the s-wave (or p-wave 
for fermions) component of the interaction, but also to 
higher partial waves components. The higher partial 
wave contributions become increasingly important as 
the scattering energy increases. To address this point, 
Giannakeas et.ai [33] recently developed a method 
involving all partial waves, and pointed out a correction 
to the position of the CIR that derives from the d-wave 
component. 

Again, extensive theoretical work has been carried 
out for systems with 2D isotropic harmonic confinement 
in the transverse plane, while other transverse confin- 
ing potentials have hardly been explored. The aim of 
the present work is to develop a systematic description 
of quasi- ID scattering in arbitrary transverse confine- 
ment. We apply the local frame transformation method 



[33], and solve the divergence problem occurring in this 
method. This divergence also arises in many other ap- 
plications of this method, analytically and numerically, 
but we believe that the systematic discussion presented 
here yields valuable insights. To illustrate our implemen- 
tation of the method, we treat the CIR in a square well 
transversely-confined system as an example. The devel- 
opment of our regularization method is inspired by the 
s-wave regularized delta function, and can be generalized 
to higher partial waves interacting through short-range 
potentials. 



length in the case of harmonic confinement, and the box 
side length for the square well confinement). H±_ is the 
"transverse" part of the Hamiltonian, and its eigenval- 
ues define the "channel threshold energies". Vj n f is the 
short range (r e is its range) spherically symmetric po- 
tential, which is not separable in Cartesian coordinates 
for which the boundary conditions are separable. In the 
remainder of this article, unless otherwise stated, a sys- 
tem of units is adopted for which h = 1 , and the reduced 
mass is /i = 1. 



II. REGULARIZED FRAME 
TRANSFORMATION 

A physical system can often be divided into different 
regions, where the controlling physics varies. In each re- 
gion, because of the difference in the symmetry of the in- 
teractions or the boundary conditions, the good quantum 
numbers can differ from one region to another. Specifi- 
cally, the quasi-lD system can be divided into two region: 
r <C L± and \z\ ^> L±, and the characteristic length L± 
in the transverse confinement is the relevant length scale 
to divide this system. At large distance in any particular 
open channel (\z\ 3> L±), the outgoing wave vector k z 
is an asymptotically good quantum number. Thus the 
phase shift for each k z asymptotically approaches to a 
constant. However, at very short range near the center of 
the confinement, where a scattering event occurs, the in- 
teraction potential has the full 3D spherical symmetry. In 
order to connect the short range scattering event and the 
large distance observables, we need to project the scatter- 
ing information near the spherically symmetric scatter- 
ing center onto the transverse channels to yield quasi- ID 
scattering information in the z-direction. The determina- 
tion of this projection is the main goal of the local frame 
transformation method. The term "local" in this con- 
text was originally coined by Fano [33] to distinguish this 
method from the usual unitary class of frame transforma- 
tions in [35] |3S] . Later, this method was extended and 
applied to several scattering problems 34, 37 39 , and 
in many cases it has resulted in excellent agreement with 
experimental observations and higher efficiency than pre- 
vious theories. Note that whereas Zhao et al. [30] claim 
to have identified a flaw in the Fano technique for trans- 
formation of the solution that is irregular at the origin, 
the present application is tested numerically in Sec. Ill 
below and it exhibits no sign of any such difficulty. 

For the Hamiltonian considered in this problem, 



H = H ± (x,y) 



2/i dz 
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the frame transformation method can be applied, be- 
cause the Hamiltonian has two clear limits: r -C Lj_ 
and \z\ 3> L±, where L± is the characteristic length in 
the transverse (confinement) plane (it is the oscillator 



This model Hamiltonian is exact for the problem of 
one light particle that scatters from an infinitely mas- 
sive particle located at the trap center. When it comes 
to treating two interacting particles (short range interac- 
tion), this treatment might still be a reasonable approx- 
imation, or at least a first step in a more general theory 
that also accounts for nonseparability of the relative and 
centcr-of-mass (COM) degrees of freedom. 



The eigenfunctions at short distance r c « r « L± can 
be represented as linear combinations of non-interacting 
channel functions: 



* 



hne 



( r ) - 2^ Fi> m i e 6iji5 mjm i - Gii m i e K Vm , lm {e), (2) 



I'm' 



where the reaction matrix K 3 encapsulates all the scat- 
tering information caused by the potential, I is the orbital 
angular momentum quantum number, m is the orbital 
magnetic quantum number. F and G are the energy- 
normalized regular and irregular solutions in spherical 
coordinates for a given scattering energy e: 



JWr) = Y,„(i)Mr) 1 = (/— j,(kr) 

s < 3 > 



Gzme(r) = Y lm (v)g le (r) g 



ni(kr), 



/ and g are the energy-normalized spherical Bessel func- 
tions, / is regular at r — ¥ 0, g is irregular at r — > 0, 
and they oscillate as r — > oo with a 7r/2 difference in 
phase. The ^-parity quantum number of the solutions is 
(— 1) , which is important to keep in mind, since later 
the 3D solutions will be projected onto outgoing waves in 
the z-dircction. The 2D isotropic harmonic confinement 
is a special case, in which the m quantum number is con- 
served in both the spherical and cylindrical symmetries, 
so there is no need to explicitly specify m [26j . However, 
for a 2D transverse trap in general, all I— and m— values 
could be coupled. In the unitary limit, of course, with 
short-range interactions, the lowest one or two partial 
waves in each symmetry are expected to dominate the 
scattering observables. 



The regular ip and irregular \ eigenfunctions at large 



distances L± <C \z\ are: 

ll>n m ,n y ,q( V ) = ^n x ,n v (x, y)(nq)~ 5 

Xn,,™ y:9 (r) = 0„ x ,„ B (x,y)(7rg)-5 



cos(qz), 
sin(^), 

sin(q\z\), 
-sign(z) cos(gz), 



(4) 

where n x , n y are quantum numbers in the trans- 
verse plane. The eigenstates in the transverse plane 
for the even parity of immediate interest here are 
cos( ^ "2 ^ a:)cos( "^ v)- The confinement in the 
x — y plane forces the scattered waves to be asymptoti- 
cally directed along either the positive or negative z-axis. 



The expression of the wave function near the scatter- 
ing center (r c < r <C L±) onto the asymptotic wave 
function (\z\ S> L±) within a fixed energy shell e can be 
accomplished by projecting the eigenfunctions near the 
scattering center onto the eigenfunctions in the asymp- 
totic region: 

Vn x ,ny,cin l: ,,ny\ r ) = 2~il,m ^n m ,n y ,q rl!ciny ;l,m,ei'lme\J c )t 

(r) = V (U- 1 ) T 



"X.n x ,ny,q nst , iri 



,q n „ ■,l,m,e* J l™zv r ) , > 

' ^ (5) 

U is the transformation matrix between the two sets of 
eigenstates, and the U_ matrix is energy dependent, since 
the expansion is performed on the energy shell. The sum- 
mation includes only solutions of the same z-parity. In 
a 2D square well trap, for which we consider only the 
s-wave interaction from the infinite-mass particle in its 
center, the projection matrix element from the transverse 
channels to the spherical s-wave is a straightforward gen- 
eralization of the spherical-to-cylindrical frame transfor- 
mation element derived in [3S], namely: 



U, 
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(6) 
where q nx .„ is the scattering momentum in each channel 



ti'k 2 



h ii 



u — , — 2u' ny ~^E ni ,n 2 ; an d we choose the branch of 
(—1)3 so that it is exp(™), 51 denotes the spacial angle in 
spherical coordinates. A similar expression arises for the 
case of 2D harmonic oscillator confinement for zero-range 
s-wave interaction: 



^n,m— 



0,q;l=0,m=0,e 



J dn^ m=0q (r)Y 00 (r) 



j (kr)^/2k/^ 
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ka± \2wj_ 
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in which ip H is the regular solution in the 2D harmonic 
confinement. As was shown in previous papers |34 [ 137 ) . 
the local frame transformation is non-unitary. This is be- 
cause the two sets of eigenstates solve two different par- 
tial differential equations (actually, in the present case, 
they solve the PDE with the same Hamiltonian but differ- 
ent boundary conditions), which are only approximately 
equal to each other in the region r c < r <^ L± |34| . 
Subsequently a number of studies [Ml EH ESS IUJ have 
discussed the resulting transformation of the reaction 
matrix K, namely: K lD = {K nin2 . ni , n2 ,} = U T K 3D U. 
The physical K 1D matrix is deduced from closed-channel 
elimination |39j . which enforces exponential decay in all 
energetically closed channels at \z\ — » oo: K} ® phys = 
K^ + iK^{l-iK^)- l K^, in which c denotes closed 
channels, o denotes open channels. All the information 
from closed channels are encapsulated into K 1 ' ys . 
When one channel is open for each overall symmetry 
(even or odd z-parity), the poles of the physical reac- 
tion matrix K} ' ys predict the positions of each CIR, 
which are values of the scattering length and/or trapping 
geometry for which Det[l — iKcc] ~ at vanishing inci- 
dent kinetic energy in the lowest channel. For the case 
of only one open channel, the lone physical K matrix 
element is the tangent of the quasi-lD scattering phase 
shift. This quasi-lD scattering phase shift in the single 
open channel recovers the pure-ID feature in the low en- 
ergy limit (since there is no concept of "closed channels" 

in a pure- ID system), having the form tan<5 = . 

k 
Clark pointed out this fact in J32] , which is equivalent to 

an observation that a "sudden" 7r/2 change occurs in the 
low energy scattering phaseshift when just an infinites- 
imal gio is turned on. We illustrate this point in Fig. 

m 

In addition to all the general properties mentioned 
above concerning quasi- ID scattering, the short-range 
potential can in some case allow further simplifica- 
tion of the problem. If K 3 is block-diagonal and 
only one angular momentum partial wave dominates 
the 3D scattering phase shift, a good approxima- 
tion to the K matrix is Kf® llm , = Sij'5 m , m > tan(<5/), 



lD,Phys 
-ID 



(E) 



K£[l - iX c (E)}-\ where A c 



TrK_ cc = tand; ^closed channels (^™*,n H ,g>i x ,„ y ;l,m,e) 
[26) . Under this single partial wave approximation, each 
pole of the physical K} ' ys is determined by the van- 
ishing of 1 — i\ c (E). 

In the treatments of this quantity X C (E) in [26l|33], its 
epxression takes one of the following forms as an infinite 
summation: 



A i _1 = Jr T,n=i \l n+ 2~ e (P" wave ) 



w=o 



«x 



n=l 



(s-wave) . 



(8) 



Moreover, in the s-wave case, Olshanii [35] encountered 
this same summation of squared transformation matrix 
elements, which is clearly divergent. A regularization 
method is developed in the following paragraphs, moti- 
vated by other analyses that have regularized different 
types of zeta functions in various contexts. 

Observe that in the definition of the Hurwitz zeta func- 
tion: 



(h(s,q) = Y^ 



1 



=*(«• 



n)' 



:,(*>!) 



(9) 



the same functional form is present in the region s € 
(l,oo) where the series is convergent. So the analytical 
continuation of the Hurwitz zeta function from s € (1, 00) 
to s G (—00, 1] can potentially regularize the divergent 
summation and yield a physically relevant value. A sim- 
ilar spirit has been implemented in the calculation of the 
Casmir force between two infinitely large planes and be- 
tween other shapes of conductors [UJ SI] • The extensive 
work on the Hurwitz zeta function in the mathematics 
community suggests a way to develop the regularization 
in general. For the present quasi-lD scattering system, 
we propose a regularization procedure that has been uti- 
lized in the mathematics community. The summation 
can be viewed as the following limit process: 



Atf&e) 




n=0 



exp(— £\/n + 1 + e) 
y/n+l + e 



(10) 



This summation is not uniformly convergent, which 
means that interchanging the order of taking the limit 
£ — > and summing over all terms might not be justi- 
fied. For any finite £, the convergence of the summation 
is guaranteed by the exponential suppression at large n. 
For £ = 0, however, the summation assumes the same 
form encountered in the I = case. The divergence in 
the local frame transformation approach can be viewed 
as having originated in an unjustifiable interchange of the 
order of these two operations. 

Olshanii et.al [2]] [S3] encountered the same divergence 
problem while solving for the scattering amplitude in the 
open channel. They treated the divergence as follows. 
The summation in Eq. [lO] can be expanded as a power 
series in £ (where £ = ^-) near the energy of the lowest 
threshold, giving: 



A H (Z,e) = F(t;)+A H (Z,e). 



(11) 



Here F is the integral approximation of the summa- 
tion right at threshold, F(£) = / °° dv c ^ ( ^/ T/) = § = 

X/s=i /"-I rffr' 6 /~~ ' wmcn exhibits the divergence 
with respect to £. A# is the regular, physically rele- 
vant part of Ah- At an arbitrary scattering energy, the 
integration approximation becomes f, dv- 



fU+l 



tion in EqJlO] can be cast as: 



A H {ti,e) = -+A H (Z,e) 



2 2 exp(-gyT+l) 



lim A r^ 00 (--(cxp(-<eVl + e) - cxp(-£y/N + e))+ 



w exp(-£Vn + e) 

2^n=l / „ 1 , ) 



= £-i00 

in which: 



+ Ah{£, — > + , e) + terms that vanish at £ = 0, 

(12) 



- 1 
A H (£ -> 0+ e) = km V -== - 2y/W+l. (13) 

n— 1 v 

The present study applies a more general concept to 
the divergent summation that has been developed in the 
mathematics community. Our treatment does not rely on 
the s-wave nature of the scattering, and can be directly 
deduced from the frame transformation result. The fol- 
lowing illustration of our method is developed specifically 
for the square well transverse trap geometry, because 
there are two quantum numbers in a general separable 
2D confinement, and the 2D square well trap is a more 
general example with quite different characteristics than 
geometries for which this type of analysis has previously 
been worked out. In addition, we show that the isotropic 
2D harmonic oscillator geometry (oj x — uj y — u) p ) is a 
special case in our treatment. 

The summation of squared [/-matrix elements in Eq. 

El 



(«x,n H )/(0,0) 

E 



1 



"■■"■' n ' ! {n x + -) 2 + (n y + -) 2 + e 



(14) 



which must be evaluated before differentiation, takes the 
following form: 



A B(4,£J - l^n^n v >0 



<"M>(-;.\ {n x + -) 2 + {n y + -) 2 + e) 



K + ^) 2 + K + 2 )2 + e 



-Y 



cxp(-Cy (n x + -) 2 + (n y + -) 2 + e) exp (_^i + ( 



{n x + -) 2 + {n y + -Y + e 



+ € 



(15) 

^(exp(— £y/l + e) — exp(— £y/N + e)). Thus the summa- The summation is over all closed-channel quantum num- 



bers in one quarter of the 2D (n x ,n y ) plane. 



The number of points in one ring \{2n — l) 2 < |n + 
(|, i)| 2 < \{1n + l) 2 can be understood as the density 
of states in the quantum number space, n = (n x ,n v ). 
The summation over all the quantum numbers can be 
approximated by an integral over the density of states in 
the 2D quantum number space, and this approximation 
can be used to separate out the singularity nature of the 
infinite sum. The integral can be evaluated by changing 
to the polar plane: 



f f dn x dn„, 
j — oo j — oo x y 



exp(-Zy/(n x + \) 2 + (n„ + \) 2 + e) 
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2irrdr 



exp(— £-\A* 2 + e) 
\/r 2 + e 
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FIG. 1. Plot of pink point n = (n x ,n y ) in the 2D integer 
plane. The ring region is defined by the range (2n — l) 2 /2 < 

with n = 3. The blue point is at 



|n + ( 
(- 



1 IS 
2' 2/ 

2' 2/' 



< 



(2ra+l) 2 
2 



27T 



limAr^oo — -(exp(-^| + e) - exp(-£ v , 



(gjV+l) 



+ e)) 

(16) 
By using the series expansion of A^ in terms of £, which 



has the same spirit of Eq. 12 we have: 



. ■ /* s 2-7T . 27T 2n . . r, xx 
4A B (e, e) = — + (-— + — cxp(-^i + e)) 



lim 



N-too 



V 



(2JV+1) 2 



+ C 



))■ 



• A r 



En=iEi < , n+( i lji^fia+a! 



(17) 



exp(-^/K + ±) 2 + (n y + \f + e) 



(n. + 5) 2 + K + |) 2 + e 



The terms inside the limiting process liniAr^oo can be 
rearranged as two parts, one of which is constant as £ — > 
0, while the other vanishes as £ — >• 0. 



ized summation: 



4A E (^0+,e) = E;=iE i 



K + o) 2 + K+o) 2 + ' 



<|n+(|,|) 



-2tt 



(2JV+1) 2 



(19) 
The ring region in {n x ,n y } is demonstrated in Fig. [I] 
The definition of the term A#(£ — > + , e = 0) is unique, 
because if one adds a term proportional to £, it will vanish 
as £ — > + , while the other term 2-7r/£ is the integral 
approximation to the summation at the threshold energy 
(e = \ in this case), which also has a unique definition. 
For 2D harmonic confinement, by the way of contrast, the 
summation over magnetic quantum numbers is trivially 
carried out because of the symmetry of that problem. 
And from this analysis, the position of the CIR in this 
geometry is calculated to be: 



M 



i 



L ± 4A £ (£^0,e) 



(20) 



2tt 
4A B (£, e) = — +4A E (£ -> + , e)+tcrms vanish at £ -)• 0+ 

(18) 
At this point, the form of the Ae function resembles that 
of the Ah function which arises for the harmonic oscilla- 
tor trap. The corresponding residual part is the regular- 



A s (£ — > 0, e) turns out to be a special case of the Ep- 
stein zeta function [35] (chapter 1, section 2.2, also briefly 
discussed below in the Appendix), a sc (e) is the energy 
dependent 3D scattering length. Many applications of 
the Epstein zeta function in physics are summarized in 
[46] . including its application in the zeta regularization 
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FIG. 2. The CIR position obtained from finite range inter- 
action model calculations is plotted along with the extrap- 
olation to zero range. The red dots are from the potential 

r 2 
V — Vo exp( -j), while the blue dots are from the poten- 



tial V = Vo 



1 



cosher 



The :r-axis is the effective range of 



the potential at the position of the CIR. The yellow point at 
r e ff — is the predicted CIR position from our regularized 
frame transformation method. The inset shows the pattern of 
convergence to the analytical value, which is seen in the nu- 
merical calculation near -4^- ~ 0.01 (the arrows indicate the 
better converged numerical values of the CIR position). At 
larger values of the 3D effective range, the numerical calcula- 
tion begins to deviate from the analytical prediction. This is 
because the range of interaction do in those calculations has 
become comparable to the CIR scattering length, and hence 
the "short" range approximation of the model potential has 
become less accurate. 



method of high energy physics J47J [48] . 

To test our proposed regularization method, we have 
first applied it to the Hurwitz zeta function summation. 
This test has verified that it gives the same numerical 
value for every tested value of q and s in agreement with 
other definitions [H3 I5U] of the Hurwitz zeta function, 
analytically continued into the region s e (—00, 1]. 



III. NUMERICAL CALCULATION OF THE 
LOW ENERGY ID SCATTERING PHASE SHIFT 

Our predicted position of the CIR can be tested further 
by performing a variational i?-matrix scattering calcula- 
tion for the square well 2D trap system using a set of B- 
spline basis functions. The CIR position for a zero-range 
3D interaction is extracted by taking its limiting value 
from finite range model potential calculations. Moreover, 
the leading finite range correction to the CIR position 
is also extrapolated and compared with the prediction 
of effective range theory. For short range interactions 
d <C a sc (CIR), the numerical model calculations agree 
accurately with the analytical prediction, as is shown in 



Fig. [2] The correction plotted versus the effective range 
of the potential exhibits a linear behavior, and the slope 
agrees with the prediction from a finite range expansion 
of the energy dependent scattering length. Some points 
deviate from the linear relation, as in Fig. [2j The reason 
for this is that the width do of the interaction is compa- 
rable to or larger than the CIR scattering length, which 
sets a new characteristic length scale that limits the ap- 
plicability of the "short" range model potential approxi- 



mation: 



dp 



«1. 



a BC (CIR) 

The yellow dot plotted at r e // = in Fig. [2] is 
the prediction from the present regularized frame trans- 
formation method, and it agrees with an extrapola- 
tion of the numerical calculation to a zero-range po- 
tential. The intercept, which corresponds to the res- 
onance position at zero energy, on Fig. [2] occurs at 

T 

^— = 5.850(±0.005) from extrapolation of the nu- 

a sc (0) 

merical results, and 5.864(±0.008) from present regular- 
ized analytical summation. The slope of the extrapo- 
lation line can be deduced from the low-energy effec- 
tive range expansion of the s-wave scattering length: 

L ± _ L 

a sc {E) 



= (0) 



L\ \ r -j^- , in which Lj_ is the width 



of the square well, and %- = E is the first threshold en- 
ergy of the 2D confinement, a sc (E) is the energy depen- 
dent scattering length, r e // is the effective range of the 
potential, which usually changes very slowly with a sc (0). 
In addition, the quasi-lD scattering phase shift recov- 
ers the behavior of the pure- ID scattering phase shift at 
perturbative values of the 3D scattering length. Fig. [3]is 
our numerical calculation for the square well transverse 
trapped system. The perturbative region {a^ol *C Lj_ is 
far away from the CIR region, and thus the transverse 
degree of freedom becomes effectively "frozen" . 



IV. BOUND STATE ENERGY FROM THE 
FRAME TRANSFORMATION METHOD 

As is well known from quantum defect theory pM 141 j . 
the reaction matrix K can be used to calculate bound 
state energy and wave function. A key question is 
whether the K matrix obtained from the local frame 
transformation can adequately predict the bound state 
energies in the quasi-lD system. As we show in the fol- 
lowing paragraph, the local frame transformation method 
can also be applied to such a bound state calculation. 

In a calculation of true bound states, all the channels 
are closed, whereby the wave function outside the range 
of the potential is: 



$ 



[MF^-K^), 



(21) 



1,3 



in which i denotes the transverse channel indices, j de- 
notes the separate linearly-independent solutions, c de- 
notes closed channels, F and G are regular and irregular 
solutions, AjS is the vector of linear expansion coefficients 



<c 








<<= 



0\ 






o- 
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FIG. 3. Plot of the quasi-lD scattering phase shift in the 
transverse square well trapped system, as a function of the 3D 
scattering length a^o- The effective ID interaction strength 



in this quasi- ID system is 



4a 3o/ L l 



B(e) 



L ± E l-B(±E )4a 3D /L ± 

is the regularized Epstein zeta function discussed in the last 
section. Eq = 5- -^- . The strips of curves clustered together 

are calculations for a range of energies in different regions of k: 
k ~ 0.1, 0.01, 0.001. The inset is the pure-ID scattering phase 
shift as a function of Q\d, which comes into the effective ID 
Hamiltonian as Vint = gioS(z). 



for the jth independent solution. Since a the bound state 
energy lies below every threshold, both F and G will con- 
tain exponential divergent components. After imposing 
the asymptotically-decaying boundary condition at in- 
finity, we find the equation that must be satisfied by the 



(iI + K)A = 0, 



(22) 



in which I is the identity matrix, and A is n x 1 col- 
umn vector of Aj . In order for this to yield a non-trivial 
solution to this system of homogenous equations, the de- 
terminant of the matrix (—il — K cc Ys must vanish. The 
closed-channel reaction matrix K\® can be derived from 
the local frame transformation, as is shown above. So 
in the quasi- ID system with only one nonzero scattering 
phase shift, the determinant is equal to: 



det | — il— K cc \ = —i(l~itain5i 



E 

i=all channels 



Uf ) (23) 



As an example, for the s-wave bound state, the relation 

between the bound state energy (E = E% ^~ L ) and the 

scattering length is: 

1 8 »-^ 1 p~ fe i|z| 



0z 



N? 



(24) 
in which Ni is the normalization constant of channel i 
in the a;j/-plane, corresponds to channel wave function fa 
and channel energy J5». This results in a similarly diver- 
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FIG. 4. Plot of the bound state energy of the weakly bound 
state in a quasi- ID system that has a 2D transverse harmonic 
oscillator trap in the transverse directions. The dots of differ- 
ent colors are from different 3D short range potentials. The 
solid line is from the analytical prediction of our local frame 
transformation treatment. The intercept zooms in to the re- 
gion where the "transition" point of the tilting behavior hap- 
pens. 



gent summation as we found in working out the quasi- ID 
scattering phase shift. Thus a similar regularization pro- 
cedure is applied to find the energy of a weakly bound 
state as a function of the 3D scattering length. 

The resulting prediction is compared with numerical 
calculations, shown in Fig. [4] Bound state energies agree 
with the analytical theory as well, but the finite range 
correction pattern is more complicated. When the bound 
state energy approaches to the threshold (1.0 Huj± in Fig. 
ffl, the finite-range potential gives a lower energy than 
does the pure zero-range result. When the bound state is 
relatively deep, the finite-range potential yields instead 
a higher energy instead. This "tilting" deviation gets 
weaker as the range of the interaction gets smaller. This 
can also be understood from the finite range expansion. 
The finite range correction term changes sign precisely 
at the point where the bound state energy changes sign, 
which can be observed on the plot. 



V. CONCLUSION 

A general method has been formulated in this study 
whereby the local frame transformation method can be 
applied to a quasi-lD scattering process and to corre- 
sponding bound state problems. The formulation pre- 
dicts the CIR position for a general transverse confine- 
ment. Our proposed procedure has been tested on a sys- 
tem with 2D square well confinement in the transverse 
directions, as well as for 2D oscillator transverse confine- 
ment where we agree with previous studies. The value of 
a L i \ = 5.864(±0.008) for the confinement-induced res- 
onance with a hard square-well confinement in the trans- 
verse direction is, interestingly enough, very close to a 



_ Oho 

a sc (0) 



factor of 4 larger than the value of -— 

a sc {0) 

1.4603... that was obtained by Olshanii for the 2D trans- 
verse harmonic confinement. As far as we can tell, this 
simple relationship is merely a coincidence. 

It should be possible to extend this regularization 
method to general 2D confinement systems having var- 
ious boundary conditions, and this direction is ripe for 
exploration in the future. 
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_ tion in 2D geometry is: 



(Bl) 
The energy normalized spherical wave function is: 



Flm,k( T ) = Y im{9,<t>)\ —3i(kr) 



(B2) 



The {/-matrix elements we use is the special case for m = 
0, I = 0. For the angle part of the integration, formula 
7.333 in [SU] is applied. 



{/n,m=0,g;i=0,m.=0, 

J~/ k 2 

ka± \2hiu± 



/dnKoo(f)< ro , g (r) 



j (kr)^2k/w 



(2n+i; 



(B3) 



Appendix A: Epstein zeta function 



The generalized form of the Epstein zeta function 
in the first chapter is: 



0)</> = z 



9i ■■■ 9 P 
hi ■ ■ ■ h p 



00, 



^ (Al) 
in which p is a positive integer, g and h are p- 
dimensional real vectors, gt,hi £ R,rn is a p-dimensional 
integer vector, m.; £ Z . 4>(x) is quadratic form of vector 
x, 4>{x) = Yjn v CpvX^Xy, in which c is p x p non-singlular 
symmetric matrix associated with %j>. The salar prod- 
uct of p-dimensional vectors (~g,h) = YZ=%9^"- e 
is the "inhomogeneity" of the generalized Epstein zeta 
function. The type of inhomogenous Epstein zeta func- 
tion we use in this paper is the special case when s = 1, 



0, namely: 



_ l _i 

2 2 





a\ 



in which <f>(x) = Eu=i x u m our case- 



Appendix B: U matrix elements 

1. 2D harmonic confinement 

In this case, at the radius r we do the frame transfor- 
mation, the energy normalized even parity channel func- 



2. 2D square well confinement 



In this case, at the radius r where the frame transfor- 
mation is performed, the energy normalized even parity 
channel function is: 



& 



n x .n y ,q 



L 



2 (2n x + l)7ra . A2n y + l)ny i 

cost- )cos(- — — — - — )(7rg) 2 cos{qz) 



L, 



Li 



(B4) 

We use the spherical expansion of plane wave to get the 
[/-matrix elements: 



(A2) exp(zk • r) = £* m i l 4nrY^(i)Yi m (r)ji(kr), 

COs(k x x) COs(k y y) COs(qz) = -^^x+ikyy+iqz _|_ e ik x x-ik y y+iqz 

8 

\ p ik x x+ik y y—iqz _i_ p ik x x—ik y y—iqz r p — ik x x+ik y y-\-iqz \ 
p —ik x x—ik y y-\-iqz _i_ p —ik x x—ik y y—iqz i p —ik x x-\-ik y y—iqz\ 



The energy normalization delta function S(E-E') comes 
into the following integration as well, so the on-shell U- 
matrix elements can be obtained by deviding the delta 



function too. 



U, 



n x ,n y ,g;/— 0,m— 0,e 



f <KlYoo(r)il> nx>nytq (r) 



j (kr)^/2k/^ 



1 Hi S d£lYoo(r)(r)e 



ik.j r 




- K + ^f - K + 2? 



(B6) 
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